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It is believed that soon after the Planck era, space time should have a semi-classical 
nature. According to this, the escape from General Relativity theory is unavoidable. 
Two geometric counter-terms are needed to regularize the divergences which come 
from the expected value. These counter-terms are responsible for a higher derivative 
metric gravitation. Starobinsky idea was that these higher derivatives could mimic 
a cosmological constant. In this work it is considered numerical solutions for general 
Bianchi I anisotropic space-times in this higher derivative theory. The approach is 
"experimental" in the sense that there is no attempt to an analytical investigation of 
the results. It is shown that for zero cosmological constant A = 0, there are sets of 
initial conditions which form basins of attraction that asymptote Minkowski space. 
The complement of this set of initial conditions form basins which are attracted 
to some singular solutions. It is also shown, for a cosmological constant A > 
that there are basins of attraction to a specific de Sitter solution. This result is 
consistent with Starobinsky's initial idea. The complement of this set also forms 
basins that are attracted to some type of singular solution. Because the singularity 
is characterized by curvature scalars, it must be stressed that the basin structure 
obtained is a topological invariant, i.e., coordinate independent. 
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I. INTRODUCTION 

The semi-classical theory consider the back reaction of quantum fields in a classical ge- 
cmetHc background. It beg., about forty years ago „.th De Wttt bi. and since then, .ts 

consequences and applications are still under research, see for example 

Different of the usual Einstein-Hilbert action, the predicted gravitational action allows 
differential equations with fourth order derivatives, which is called the full theory [l], see 
also |3|. 

There are several problems connected to the full theory which prove that it is not consis- 
tent with our expectation of the present day physical world, see jJI and references therein. 
In this work, Simon and Parker-Simon suggest a perturbative approach which is very inter- 
esting. 

n 

The full higher order theory was previously studied by Starobinsky j^], and more recently, 
also by Shapiro, Pelinson and others j^. Starobinsky idea was that the higher order terms 
could mimic a cosmological constant. In only the homogeneous and isotropic space time 
is studied. 

The full theory with four time derivatives is addressed, which apparently was first inves- 
tigated in Tomita. article R for general BrancM I spaces. They found that the presence 

of anisotropy contributes to the formation of the singularity. Berkin's work shows that a 
quadratic Weyl theory is less stable than a quadratic Riemann scalar E? . Barrow and Hervik 
found exact and analytic solutions for anisotropic quadratic gravity with A > that do not 
approach a de Sitter space time. In that article Barrow and Hervik discuss Bianchi types II 
and Vlh and also a very interesting stability criterion concerning small anisotropics. There 
is also a recent article by Clifton and Barrow in which Kasner type solutions are addressed. 
H. J. Schmidt does a recent and very interesting review of higher order gravity theories in 
connection to cosmology |8|. 

On the other hand, numerical investigations of initial conditions, was already addressed 
in Einstein's theory of gravitation for the homogeneous Bianchi IX solutions. Basin of 
attraction to asymptotic solutions, and the fractal nature of the basin boundary are strong 
indications that the early stages of the Universe were very chaotic, according to Einstein's 
theory, see for example 9|. 

For general anisotropic Bianchi I homogeneous space times, the full theory reduces to a 
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system of nonlinear ordinary differential equations. The numerical solutions for this system 
with derivatives of fourth order in time was previously obtained by us lol |. 

In this present work a numerical analysis of initial conditions for this nonlinear system 
mentioned above, is given. 

Only the vacuum energy momentum classical source is considered in the full theory. It 
should be valid soon after the Planck era in which vacuum classical source seems the most 
natural condition. The chaotic nature of the singularity is not addressed. 

There is not a unique way to totally characterize a singularity. In this numerical investi- 
gation, the singularity is defined when RabR"^^ > 10^. 

It is found that for zero cosmological constant A = 0, there are basins that are attracted 
to zero, constant 4-curvature solutions. This proves that according to this theory, there are 
regions in the space of initial conditions, for which the Minkowski geometry is structurally 
stable. That is, for considerably anisotropic initial conditions isotropisation occurs, and 
Minkowski space is obtained asymptotically, without the need for fine tuning. 

On the other hand for a nonzero cosmological constant A > the basins attract to a 
specific de Sitter geometry, which is a constant homogeneous, non null 4-curvature solution. 
Thus indicating that Starobinsky's inflation is structurally stable in the above sense. 

Also, the complement of the set of initial conditions for the stable, physically accepted 
solutions, for both cases A = and A > 0, form basins of attraction to singular, in some 
sense. Universes. Thus physically acceptable initial conditions evolve in very few Planck 
times, to Universes with very high curvature scalars. In this sense, this theory certainly is 
not a complete one. 

It should be mentioned that the isotropisation process for non zero cosmological constant 
is not a peculiarity of these higher derivatives theories. It occurs in ordinary Einstein 



General Relativity theory. This is proved in the very interesting article [ll| , where for A > 
all Bianchi models, except for highly positively curved Bianchi IX become asymptotically 
de Sitter. The remarkable difference between General Relativity theory and this higher 
derivatives theories is that the isotropisation depends on the initial conditions. Isotropisation 
is strongly dependent on initial conditions for quadratic curvature theories such as (P). 

The following conventions and unit choice are taken i?^^^ = P^^^ — Rat = Racbi ^ai 
metric signature — h ++, Latin symbols run from — 3, Greek symbols run from 1 — 3 and 
G = h = c=l. 
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II. THE MODEL 



The Lagrange function is, 



Metric variations in the above action resuhs in 



A + R + a[ W^Rab -^R^]+ (3R^ 



+ (1) 



Eab = Gab + - H^J + aH^ab - Tab - \l^gab. (2) 



where 



H^b = \9abR^ - ^RRab + 2R.ab ' ^DRgab, (3) 

-^afe'* — lj^9abRmnR'^"' + R;ab — '^R'^^ Rcbna — D-Rafe — l^Rdabi (4) 
Gab = Rab — -QabR, (5) 

and Tab is the energy momentum source, which comes from the classical part of the La- 
grangian Cc- Only vacuum solutions Tab = Cc = will be considered in this paper since it 
seems the most natural condition soon after the Planck era. 

The covariant divergence of the above tensors are identically zero due their variational 
definition. The following Bianchi Type I line element is considered 

ds'^ = -dt^ + [ai{t)fdx^ + [a2{t)fdy'^ + [a^{t)fdz^, (6) 

which is a general spatially fiat and anisotropic space, with proper time t. With this line 
element all the tensors which enter the expressions are diagonal. The substitution of © in 
((2|) with Tab = 0, results for the spatial part of Q, in differential equations of the type 



d"" 


= /i 


f d^ 




ydt^ 




= f2 


} d^ 




ydt^ 


d^ 


= h 


} d^ 


dF""' 


,dt^ 



ai, di, di, ai ) (7) 

(8) 
(9) 



where the functions fi involve the ai, 02, 03, and their derivatives in a polynomial fashion, 
see the Appendix El The very interesting article [l2] shows that the theory which follows 
from (0) has a well posed initial value problem. This question has some similarities which 
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General Relativity theory since it is necessary to solve the pro blem of boundary conditions 



12 1, once a consistent boundarA 



and initial conditions simultaneously. According to Noakes 
condition is chosen initially, the time evolution of the system is uniquely specified. In |l2| 
the differential equations for the metric are written in a form suitable for the application of 
the theorem of Leray The numerical solutions of the partial differential equations in 
(121) are not going to be addressed in this present work. 

For homogeneous spaces the differential equations in (j21) reduce to non linear ordinary 
differential equations. Then, instead of going through the general construction given in 
jl^ . in this particular case, the existence and uniqueness of the solutions of ((H)-© reduce 
to the well known problem of existence and uniqueness of solutions of ordinary differential 
equations. For a proof on local existence and uniqueness of solutions of differential equations 



see, for instance 

Besides the equations ©-Q, we have the temporal part of (j21). To understand the role 
of this equation we have first to study the covariant divergence of the equation Q, 

Remind that the coordinates being used are © and that Tab = 0. Since the differential 
equations (jZI)-© are solved numerically E""^ = 0, 

doE'' + T:,E'' + T',,E'' = 0, 

E^\t) = 0, (10) 

where Eqq is the 00 component of 0. If Eqq = initially, it will remain zero at any instant. 
Therefore the equation Eqq acts as a constraint on the initial conditions and we use it to 
test the accuracy of our results. 

For a space-like vector f " and a time-like vector t"" = (1, 0, 0, 0) tidal forces are given by 
the geodesic deviation equations 



-^00/3 ~ ^cel3 



00/3' 



The theory predicted by is believed to be correct if the tidal forces are less than 1 in 
Plank units units, 

l-Roaol — -'- summation.) (11) 
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When this condition is not satisfied, quantum effects could introduce further modifications 
into (dJ. 

III. ANALYSIS OF INITIAL CONDITIONS 

Ordinary differential equations are deterministic in the sense that initial conditions specify 
the dynamical evolution of the system uniquely. 

The interest in the numerical investigation is in the sensibility upon initial conditions. In 
this work, the system is treated as an exit system. 

It is a relatively simple approach in which asymptotic classes of solutions are recognized. 
The initial conditions are identified with respect to which asymptotic solution the system 
evolves to For this purpose, holes are cut in the phase space and if the system falls into 
a specific hole, this specifies the asymptotic class. The closure of the set of initial conditions 
which asymptotes a given solution is called the basin of attraction. 

The presence of basins of attraction is an indication of structural stability, since for many 
initial conditions, without the need for fine tuning, the same class of asymptotic solution 
is obtained. It is also a topological invariant, since a smooth coordinate transformation is 
impossible to modify the Hausdorff dimension of a set jl^. 

In subsections IIII Al and ITTTB| the differential equations ((Zl)-® are numerically solved 
and the initial conditions are chosen in the following way. All the higher derivatives of the 
metric set initially to zero 

-^aa{t = 0) = 'dait = 0) = and ai{t = 0) = a2(t = 0) = as^t = 0) = 1 (12) 
dt'^ 

The Hubble constants in each direction. Appendix 1X1 

H^ = ^^0, H, = ^^,0 (13) 
together with the Hamiltonian constrain ()10|) fix the initial value of 

^3 = ^. (14) 

^3 

This value is not unique since eq. (jlUj) in the coordinates (jH)) form a polynomial of degree 
3 in the first derivatives. For arbitrary initial values for Hi and H2 there could be at most 
3 real distinct initial values for H3. One of these values is chosen. It can be seen that the 
initial values are consistent with the condition |TT|). resulting in null tidal forces. 
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FIG. 1: Numerical solution of 0-®! with the initial condition described in eqs. (|12() - (|14|) . with 
a = 1, P = —1 and A = 0, plotted against time, in Planck units. It can be seen the isotropisation 
towards Minkowski space. 
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FIG. 2: Numerical solution of 0-®, with the initial condition described in eqs. (fT2|) - (fTl|) . with 
a = 1, /3 = —1 and A = 0, plotted against time, in Planck units. It can be seen the formation of a 
singularity with Hi < and H2 < and < 0. 

A. Stability of Minkowski space 



In this subsection the cosmological constant is set to zero A = 0. The numerical solutions 
of ©-© are identified according to their asymptotic classes. In FIG. ^ it is shown the 
isotropisation process toward Minkowski space. A singularity for which all the 3 dimensions 
shrink to zero, with Hi < 0, H2 < and H3 < when t — 00 is shown in FIG. |2l In FIG. El 
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FIG. 3: Numerical solution of 0-®! with the initial condition described in eqs. (|12() - (|14|) . with 
a = 1, /3 = —1 and A = 0, plotted against time, in Planck units. It can be seen the formation of a 
singularity with at least i/i > or H2 > or > 0, and at least one of the H is negative. 

it is shown a singularity for which one of the dimensions increases, either ifi > or > 
or H3 > and at least one of the H is negative, when t —>■ 00. In FIG. @J a black, white or 
gray point is addressed to a initial condition whether it falls into the classes shown in FIG. 
d FIG. El or FIG. El respectively. 



B. Stability of de Sitter space 

In this subsection a non zero cosmological constant is chosen A = 0.02. The map of 
initial conditions for the differential equations ©-(jHl), is almost identical to the last section, 
except for the black point. Each black point in FIG. El corresponds to initial conditions 
whose solution asymptote a constant 4 curvature space R = Rabd"'^ = —0.04 shown in FIG. 
ini In FIG. inithe time evolution of a black point in FIG. El shows the isotropisation process 
toward a de Sitter, constant curvature solution. The basins of attraction to this same de 
Sitter solution, is an indication of stability in some sense. 
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FIG. 4: A two dimensional section defined by (fT^ - (fT^ and (fTU]) . for the phase space of 0-®, 
with 400 X 400 initial values for Hi and H2. a = 1, (3 = —1 and A = 0. The black regions are 
the basins of attraction to the Minkowski space. The white and gray regions are initial conditions 
which asymptote the singularities shown in FIG. and FIG. |31 respectively. 
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FIG. 5: A two dimensional section defined by (fT2|) - (|14|) and (|T?1|) . for the phase space of 0-®, 
with 400 X 400 initial values for Hi and H2. a = 1, P = —1 and A = 0.02. The black regions 
are the basins of attraction to the de Sitter space sown in FIG. IHl The white and gray regions are 
initial conditions which asymptote the singularities shown in FIG. Eland FIG. El respectively. 
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FIG. 6: Numerical solution of Q-©, with the initial condition described in eqs. (fT^ - (fTl|) . with 
a = 1, P = —1 and A = 0.02, plotted against time, in Planck units. It can be seen the isotropisation 
toward a de Sitter space with R = Rabd""^ = —0.04. 

IV. CONCLUSIONS 

In the present work it is considered general anisotropic Bianchi I homogeneous space 
times given by the line element ©. For this line element, the theory given in (^, reduces to 
a system of ordinary nonlinear differential equations with four time derivatives. This system 
is numerically investigated. In particular, an analysis of initial conditions is given, which is 
the interest from the dynamical system point of view. The supposition of a homogeneous 
Universe is artificial, and still presents a first generalization for its very primordial stages. 

Anyway, there is a well known conjecture that dissipative processes can take place in a 
infinite dynamical system reducing the number of degrees of freedom to just a few. In some 
sense, this argument reduces the artificiality of the supposition of a homogeneous Universe. 

Only the vacuum energy momentum classical source is considered in the full theory. It 
should be valid soon after the Planck era in which vacuum classical source seems the most 
natural condition. The chaotic nature of the singularity is not addressed. 

There is not a unique way to totally characterize a singularity. In this numerical investi- 
gation, the singularity is defined when RabR"^^ > 10^. 

It is found that for zero cosmological constant A = 0, there are basins that are attracted 
to zero, constant 4-curvature solutions. This proves that according to this theory, there are 
regions in the space of initial conditions, for which the Minkowski geometry is structurally 
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stable. That is, for considerably anisotropic initial conditions isotropisation occurs, and 
Minkowski space is obtained asymptotically, without the need for fine tuning. The fact that 
the Minkowski space is an attractor to general Bianchi I spaces must be connected to non 
linear effects, see Appendix 1X1 

On the other hand for a nonzero cosmological constant A > the basins attract to a 
specific de Sitter solution. That is, constant homogeneous, non null 4-curvature solutions. 
Thus reproducing Starobinsky's idea. 

Also, the complement of the set of initial conditions for the stable, physically accepted 
solutions, for both cases A = and A > 0, form basins of attraction to singular, in some 
sense. Universes. Thus physically acceptable initial conditions evolve in very few Planck 
times, to Universes with very high curvature scalars. In this sense, this theory certainly is 
not a complete one. 

It should be mentioned that the isotropisation process for non zero cosmological constant 
is not a peculiarity of these higher derivatives theories. It occurs in ordinary Einstein 
General Relativity theory. This is proved in the very interesting article , where for A > 
all Bianchi models, except for highly positively curved Bianchi IX become asymptotically 
de Sitter. The remarkable difference between General Relativity theory and this higher 
derivatives theories is that the isotropisation depends on the initial conditions. Isotropisation 
is strongly dependent on initial conditions for quadratic curvature theories_such as (P). 

The analytical solutions found by Barrow-Hervik and Clifton-Barrow [7] can be under- 
stood as limit sets in the space of solutions of the quadratic theory (^. In particular, 
Barrow-Hervik found solutions for which isotropisation does not occur. In this present work 
it is shown that the initial conditions for which isotropistation does not occur, form a set 
with non zero measure. Every gray point in FIG. 0] and FIG. El corresponds to solutions for 
which isotropisation does not occur. A bifurcation analysis with respect to the parameters 
space, a, /3, A is beyond the scope of the present work. 
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APPENDIX A: THE DIFFERENTIAL EQUATIONS 

The following line element is chosen 

With this choice, the first time derivatives of the functions Wi{t) are related to the functions 
ai{t) of (jni) in the following manner 

yi = wi = - (Al) 

ai 

0,2 , . » 

y2 = W2 = — (A2) 

Vs = W3 = — (A3) 

as 

i/4=«ii=^-f^y (A4) 

i/5=^.=^-f^y (A5) 

(A6) 

02 ' 





d's 








Z/6 = 


W3 = — 










as 












1 




a'l 


1/7 = 




ai 




ai 






1 




a2 


1/8 = 




a2 




a2 






1 




a's 


2/9 = 




as 




as 



3^-2(-lj I (A8) 

as ■ 



Sri _ 2 ( ) ) . (A9) 
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The differential equations for the quadratic theory with the classical source Tab = 0, Q, are 
equivalent to 

-^wi = -^y^ym -y 2/3^21/1 +y 1/22/12/3 + 5/3 yiysys - 13/3 2/12/22/4 

14 

-13/3 2/12/32/4 + 1/9 2/12/62/2 + ^/Syiy^ye + 2/3 2/22/32/6 - y 2/22/32/4 + 2/3 2/22/32/5 

1 2 1 2 1 2 1 

+ 1/9 2/12/32/5 ~ ( ?/32/i2/2 + Y^ 2/32/22/1 + y^^ 2/22/i2/3 + ^ 2/i2/62/2 

1 1 1 1 1 / 2 2 

+ ^ 2/22/32/4 + ^ 2/12/32/5 - Y^ 2/52/4 - Y^ 2/62/4 - Y^ 2/62/5 + 1/27 2/1 2/2 

+ 1/272/1^2/3' + 1/272/2^2/3' - 1/362/3=^2/1 - 1/362/1^2/2 - 1/362/1^2/3 - 1/36 2/2^2/1 
-1/36 2/3^2/2 - 1/36 2/2^2/3 - ^ 2/12/7 - ^ 2/22/8 - ^ 2/32/9 + y^ 2/82/1 + y^^ 2/22/7 
1 1 1 1 2 1 2 1 2 1 4 



2/92/1 + Y5^ 2/92/2 + Y5^ 2/32/8 + Y5^ 2/4 +Y5^2/5 +Y5^2/6 +Y5^2/2 

I4I2 I2 l2|l 1I4I <3-l 

2/3 - 77 2/1 2/4 - ^ 2/2 2/5 - 77 ?/3 2/6 + T7T77 2/32/7 + tttt 2/i «P 



108 54" " 54---- 54---- 108 108 

- (-1/18 2/6 - 1/18 2/5 + 1/6 A - 1/18 2/3' - 1/36 2/32/2 - 1/18 2/2' - 1/18 2/4 - 1/36 2/32/1 

11 11 5 5 

-1/18 2/1' - 1/36 2/22/1) /(C/?) - y 2/52/4 - y 2/62/4 + 1/9 2/62/5 - y^ 2/1^2/2' - —yi^ys"^ 

7 22 

+ Y^ 2/2^2/3' + 2/3 2/3^2/1 - 2/3 2/1^2/2 - 2/3 2/1^2/3 + 2/3 2/2^2/1 - y 2/i2/7 - 1/9 2/22/8 

-1/9 2/32/9 - 1/9 2/82/1 - y 2/22/7 - l/9 2/92/i + 2/9 2/92/2 + 2/9 2/32/8 - y^ 2/4' 

^2 ^2 ^4 ^4 , 2 ^2 ^2 16 

"18 ~ 18 ~ 18 ~ 18 - 1/9 2/1 2/4 - - 2/2 2/5 - g 2/3 2/6 - y 2/32/7 

-4/3 2/52/1' - 4 / 3 2/62/i' + 1/3 2/62/2' + l/3 2/3'2/5 - ((4/3 2/32/12/2' - 8/3 2/32/22/1' 

+4/3 2/22/12/3' - 4 2/12/22/5 - 4 2/12/22/4 - 4 2/12/32/4 - 4/3 2/12/62/2 - 4 2/i2/32/6 + 8 2/22/32/6 

+8/3 2/22/32/4 + 82/22/32/5 - 4/3 2/12/32/5 - 4/3 2/52/4 - 4/3 2/62/4 + 8/3 2/62/5 - 8/3 2/i'2/2' 

-8/3 2/i'2/3' + 16/3 2/2'2/3' - 4 2/1^2/2 - 4 2/1^2/3 + 4 2/3^2/2 + 4 2/2^2/3 - 8/3 2/12/7 

+4/ 3 2/22/8 + 4/3 2/32/9 - 8/3 2/82/1 +4/3 2/22/7 - 8/3 2/92/1 +4/3 2/92/2 + 4/3 2/32/8 - 8/3 2/4' 

32 

+4/3 2/5' + 4/3 2/6' + 4/3 2/2' + 4/3 2/3' - — 2/i'2/4 + 16/3 2/2'2/5 + 16/3 2/3'2/6 
+4/3 2/32/7 - 4 2/52/1' - 4 2/62/i' + 4 2/62/2' + 4 2/3'2/5 - 8/3 2/1') + (-1/3 2/3' - 1/3 2/5 
-1/3 2/2' - 2/3 2/32/2 - 1/3 2/6 + 2/3 2/1' + 2/3 2/4 + 1/3 2/22/1 + 1/3 2/32/i) /G) a'' 
+^ 2/1' (AlO) 

= -y 2/32/12/2' + y 2/32/22/1' + '^y2yiyz - 13/32/12/22/5 + 5/32/12/22/4 + 2/32/12/32/4 

14 

+ 1/9 2/12/62/2 + 2/3 2/12/32/6 + 5/3 2/22/32/6 + l/dy2y3y4 - 13/3 2/22/32/5 - -^ymy^ 
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/I 2 1 2 1 2 1 1 

~ i 108 ^^^^^^ ^ 108 ^^^^^^ ^ 108 ^^^^^^ ^ 54 ^^^^^^ ^ 54 ^^^^^"^ 

+ 77 1/11/31/5 - ^51/4 - 7^ 1/61/4 - yeVh + 1/271/1^1/2^ + Ijl^yxvi 
54 108 108 108 

+l/27y2'z/3' - 1/36 1/3 V - 1/361/1^1/2 - 1/361/1^1/3 - 1/36 1/2V - 1/361/3^1/2 

. /o. 3 1 1 1 1 1 1 

-1/36 y2 ?/3 - ^ yiy7 - ^ ^2^8 - ^ 1/31/9 + l/Sl/l + V'lVl + ^91/1 

^'^^ + ik + M + M + ik + M + M 

l/i'l/4 - ^ l/2'l/5 - ^ l/s'l/e + ^ 1/31/7 + Vx^ a/3"' - (-1/18 i/e 
-1/18 1/5 + 1/6 A - 1/18 1/3' - 1/361/31/2 - 1/18 1/2' - 1/18 1/4 - 1/36^3^1 
-1/18 yx - 1/36 1/2I/1) /(G/3) - — y^y^ + 1/9 y6l/4 - y I/6I/5 - Vxy-i 

+^ Z/i'y3' - ^ y-iVz + 2/31/1^1/2 - 2/31/2^1/1 + 2/31/3^1/2 - 2/31/2^1/3 

22 16 
-1/9 yiyr - y yiy% - 1/9 ^3^9 - — ysVi - 1/9 1/2I/7 + 2/9 ygyi - 1/9 1/91/2 

^3^8 ~ ^ - ^ ys^ - ^ 1/6^ ^ Ts^^^ ~ Ys^^^ ~^ ^^^^^ ~ ^^^^^ 

-^1/3^1/6 + 2/9 1/31/7 + 1/3 1/61/1^ - 4/3 1/61/2^ ~ ^ ~ 4/3l/2V + 1/3 1/3^1/4 
- ((-8/3l/3l/il/2^ + 4/3 l/3l/2l/l^ + 4/3l/2l/ll/3^ - 4l/il/2l/5 - 4l/il/2l/4 + 8 I/1I/3I/4 

-4/3 yiy6y2 + 8 yiy3y6 - 4 y2y3y6 - 4/3 y2y3y4 - 4 y2y3y5 + 8/3 yiy3y5 - 4/3 y^y^ 

+8/3 1/61/4 - 4/3 1/61/5 - 8/3 1/1=^1/2' + 16/3 yi^y^^ - 8/3 1/2^1/3' + 4 i/3'l/i 

+4 i/i3y3 - 4 1/2^1/1 - 4 1/2^1/3 + 4/3 1/1I/7 - 8/3 1/2I/8 + 4/3 1/3I/9 + 4/3 ysyi - 8/3 1/2I/7 

+4/3 ygyi - 8/3 ygy2 + 4/3 ^3^8 + 4/3 ^4' - 8/3 ^5' + 4/3 yg' - 8/3 ^2^ + 4/3 ^3^ 
32 

+16/3 yi% - — y2^y5 + 16/3 y^^y^ + 4/3 1/3I/7 + 4 i/ei/i^ - 4 1/61/2=^ + 4/3 1/1^ 
-4 1/2^^4 + 4 1/3^^4) (3 + (-1/3 1/3' + 2/3 1/5 + 2/3 1/2' + 1/3 1/3I/2 - 1/3 y^ 
-1/3 yi^ - 1/3 y4 + 1/3^2^1 - 2/3 ^3^1) /G) a"^ (All) 
y y3yiy2^ + y y3y2yi^ - y y2yiy3^ + 2/3 yiysys + '^/^ymV'i 

14 

+5/3^1^3^4 - y yiy6y2 - 13/3^1^3^6 - 13/3^2^3^6 + 1/9^2^3^4 + 5/3^2^3^5 

/I 2 1 2 1 2 1 

+1/9 yiy3y5 - I ^3yiy2 + — y3y2yi + — y2yiy3 + ^ yiy6y2 
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1 1 1 1 1 1^22 

+ ^ 1/21/31/4 + ^ 1/11/31/5 - 1/51/4 - 1/61/4 " 1/61/5 + 1/27 ?/i |/2 

+l/27yi^7/3' + 1/27 y2V - l/Se^aV - l/Seyi^ys - 1/36^1^7/3 - 1/36|/2V 
-1/36 ya^ya - 1/36 ya^ya - ^ yiy? - ^ ^2^8 - ^ ^3^9 + ysl/i + ^2^7 

+1^ - va y^'y' - k ^^'^^ - k y^^' + io8 



^1/31/7 + - (-1/18 ye 



-1/18 ys + 1/6 A - 1/18 ya' - 1/36 yays - 1/18 ya' - 1/18 y4 - 1/36 yayi 

-1/18 yi^ - 1/36 ysyi) /{G(5) + 1/9 y5y4 - y y6y4 - y y&Vb + ^ yi^y2^ - ^ yi^ys^ 

y2'y3' - 2/3 ya^yi + 2/3 yi^ys - 2/3 ys^^ys + 2/3 y^^'y^ - 1/9 yiyy - 1/9 ym 

22 ^ ,^ ^ ,^ 16 16 5,52 

-- g- + 2/9 ysyi + 2/9 y2y7 - y ygyi - y ^9^2 - 1/9 y3y8 - - 

- ye^ - ^ ^2^ + ^ y3^ - ^ yi^y4 - ^ y2^y5 - ^/Qys^ye - 1/9 ysy? + i/3 ysyi^ 

-4/3 ya^ys - ((4/3yayiy2^ + 4/3yay2yi^ - 8/3y2yiy3^ + 8yiy2y5 + 8yiy2y4 

-4 yiy3y4 + 8/3 yiy6y2 - 4 yiysye - 4 y2y3y6 - 4/3 y2y3y4 - 4 y2y3y5 - 4/3 yiyays 

+8/3y5y4 - 4/3y6y4 - 4/3 yeys + 16/3yi'y2=^ - 8/3yi'ya' - 8/3y2\a' 

-4 ya^yi + 4 yi^y2 + 4 y2^yi - 4 ya^y2 + 4/3 yiyy + 4/3 y2y8 - 8/3 yayg + 4/3 ygyi 

+4/3 y2y7 + 4/3 ygyi + 4/3 y9y2 - 8/3 yayg + 4/3 y4' + 4/3 y,^ - 8/3 ye' + 4/3 y2' 

32 

-8/3 ys' + 16/3 yi2y4 + 16/3 ys'ys - — ys'ye - 8/3 yayr + 4 ysyi' - 4 ya'ys 

+4/3 y/ + 4 y2'y4 - 4 ya2y4) P + (2/3 ya' - 1/3 ys - 1/3 y2' + 1/3 yay2 + 2/3 ye 

-1/3 yi' - 1/3 y4 - 2/3y2yi + l/3yayi) /G) a-' ^ ^yi' + 1/3 y2V 

-4/3y32y4 (A12) 
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